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1 

The present invention relates to a system for analysing ECG curvature where at least 
one among a number of different parameters is isolated, which system has a input 
means connected to an ECG source, where the different parameters of a received ECG 
curvature are indicated and/or isolated and for indicating possible symptoms which 
5 relates to or are indications of certain deceases, where said deceases are known to in- 
fluence the ECG curvature. 

The present invention further relates to a method for analysing ECG curvature, which 
curvature contains a number of parameters. 

10 

The aim of the invention is to achieve a system and a method for objective, fast and 
effective indication of a number of symptoms derivable from an ECG curve which 
may be indicative of one or more diseases. 

1 5 This can be achieved with the system previously described if a first number of selected 
parameters, are combined in at least a first mathematical analysis, where the result of 
the analysis can be represented as a point in a coordinate system comprising at least 
two axes where the system can compare the actual placement in the coordinate system 
with a number of reference parameters stored in the system for indicating diseases 

20 having influence on the ECG curvature. 

Hereby, it is achieved that any symptom of a disease having an indication (influence) 
in the ECG curvature can be detected in an objective, automated and very fast way. 
The system might be used under field conditions such as in ambulances or in other 

25 situations where a fast indication of heart diseases is needed in order to help the pa- 
tient in a correct way as early as possible. The analysis that takes place in an ambu- 
lance on its way to the hospital, can by transmitting the results to the hospital, allow 
the doctor at the hospital to give feedback to the personnel in the ambulance so that 
the correct treatment of the patient may start. At the same time, the hospital can pre- 

30 pare the correct activity for the incoming patient. The system could be very important 
for ECG analyses for all non-specialists in the field if they have to analyse a ECG 
curvature. 
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2 

The scope of the invention can also be fulfilled with a method for analysing the ECG 
curvature if the method incorporates the steps of: 

- receiving ECG curvature from a source, 

5 - indicating a number of different parameters contained in the received ECG 

curvature, 

- storing the parameters in storage means, 

- selecting disease specific parameters in the storage means 

- combining selected parameters in mathematical analysing means 

10 - representing the result of the mathematical analysis as a point in a coordinate 

system, comprising at least two axes, 

- comparing the actual placement in the coordinate system with a number of ref- 
erence parameters stored in a memory, 

- indicating possible diseases having the determined influence on the ECG cur- 
15 vature. 

In this way as already described, a very effective analysis of the ECG curvature is 
achieved. 

20 Below the parameters are illustrated on a typical ECG curvature in order to illustrate 
the different curve sections isolated by the analysis, and referred to by parameters, 
intervals and complexes. 
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5 The heart generates an electrical signal ECG. The waves at the ECG-signal P,Q,R,S,T 
and U are due to depolarisation and repolarisation of the heart. 

The QRS-complex contains the Q,R and S wave, and goes from the Q-onset to S- 
offset. 

10 

The QT interval starts at the beginning of QRS (Q-onset) and ends at the T-wave on- 
set. 

The RR interval goes from one R-peak to the following. 

15 

The PR interval is from the beginning of P to the beginning of QRS (Q-onset). 
QRS duration is the width the QRS complex. 
20 The ST-segment starts at the QRS-complex offset and ends at the T-wave onset. 
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The analysing process is repeated in the system for further selected parameters in order 
to indicate further possible diseases or symptoms. Hereby, it is achieved that the sys- 
tem or the method can be repeated several times with different combinations of pa- 
rameters. In order to analyse for a high number of possible diseases, each time a result 
of an analysis is achieved, the parameters are compared with stored parameters where 
every stored parameter is an indication of a symptom of a known disease. With the 
system, a deviation of parameters from the stored data indicating symptoms of an ex- 
act disease may also be interpreted for further reference. Parameters could be con- 
trolled, and if the storage means has an enormously high number of reference parame- 
ters relating to various symptoms, testing for all symptoms and correlated diseases 
having influence on the ECG curvature may be performed. 

The system or method divides the parameters into at least two main groups, which 
groups contain parameters of symmetry, flatness and duration relating to the actual 
ECG curvature. In this way, it is achieved that the parameters are grouped in the sys- 
tem, and in every group, they can be further subdivided into a specific number of pos- 
sible parameters. Keeping the number of parameters relatively small, the analysis takes 
place in a faster way. 

The group of symmetry might contain at least the following parameters: 

5 1 Skewness evaluated from Tstart to Tend, 

52 Skewness evaluated from Tstart to Tend with Ttop as mean, 

53 Skewness evaluated in a symmetric interval, 10 % of the Tstart-Tend interval 
surrounding Ttop with Ttop as mean 

54 Skewness evaluated in a symmetric interval, 20 % of the Tstart-Tend interval 
surrounding Ttop wit Ttop as mean, 

55 Ratio of the time interval from Tstart to Ttop and the time interval from Ttop to 
Tend. 

56 Ratio of the average slope from Tstart to Ttop and from Ttop to Tend. 



The group of flatness might contain at least the following parameters: 

F 1 Kurtosis evaluated from Tstart to Tend, 

F2 . Fl normalized by the absolute Rtop-Qnadir value, 

F3 Kurtosis evaluated from Tstart to Tend with Ttop as mean, 

F4 F3 normalized by absolute Rtop-Qnadir value, 

F5 Kurtosis evaluated in a symmetric interval, 10 % of the Tstart-Tend interval sur- 
rounding Ttop with Ttop as mean, 

F6 F5 normalized by absolute Rtop-Qnadir value, 

F7 Kurtosis evaluated in a symmetric interval, 20 % of the Tstart-Tend interval sur- 
rounding Ttop with Ttop as mean, 

F8 Kurtosis normalized by the value of Rtop with Ttop as mean, 

* 

F9 Ratio of the total area under the T-wave from Tstart to Ttop and the correspond- 
ing time interval, 

F 10 F9 normalized by absolute Rtop-Qnadir value, 

Fl 1 Ratio of the total area under the T-wave from Ttop to Tend and the correspond- 
ing time interval, 

F 1 2 F 1 1 normalized by absolute Rtop-Qnadir value 

F 1 3 Ratio of the total area under the T-wave from Tstart to Tend and the correspond- 
ing time interval, 

F 14 F 1 3 normalized by absolute Rtop-Qnadir value, 

F15 Ratio of the height of Rtop and the width of the Tstart-Tend interval. 

The group of duration might contain at least the following parameters: 

QTc The Q-T interval normalized by the square root of the R-R interval according to 
Bazett's formula , 

D2 Time interval from Tstart to Tend, 



6 

D3 Time interval from Tstart to Ttop, 
D4 Time interval from Ttop to Tend. 

The groups of parameters could contain further parameters and the group may contain 
5 a number of sub groups. 

Combining parameters, the system and/or method may form the different groups. 
When combining parameters from different groups, a much better result is achieved . 
than when only using parameters from the same group. The parameters can be an ele- 
10 vation of the curve; they can be the morphology of the curve; or they could be time- 
deviations as an example of possible parameters. When combining parameters, a pre- 
cise analysis can take place because a specific combination of parameters can indicate 
a specific disease and it is possible effectively to select between ECG-signals that look 
very much alike, but which indicate different diseases. 

15 

The system and/or method can analyse the QT interval of the ECG curvature for indi- 
cating Long QT syndrome. This way, the Long QT syndrome can be indicated in an 
objective and effective manner which might occur with children right after they are 
born in order to start a treatment of Long QT syndrome as early as possible. The 
20 method can differentiate between different genotypes of the Long QT Syndrome, 
which is important for the treatment. Hereby can be achieved that the correct medical 
treatments can be started. 

The system and/or method can analyse for ST-elevation myocardial infarction by ana- 
25 lysing at least the following parameters: ST elevation, ST morphology, T wave mor- 
phology and Q wave morphology. In this way, a very effective indication for ST- 
elevation myocardial infarction is achieved at the correct activity can be stated as early 
as possible. 

* 

30 The system and/or method can analyse for Non ST-elevation myocardial infarction by 
analysing at least the following parameters: ST depression, T wave morphology and Q 
wave morphology. Non ST-elevation myocardial infarction can also be detected in a 
highly effective way and correct treatment can be stated. 



Copy provided by USPTO from the IFW Image Database on 11/18/2004 



The system and/or method can analyse for Cardiomyopathia by analysing at least the 
following parameters: P wave morphology, QRS duration, S Wave morphology, T 
wave morphology. Cardiomyopathia can in this way be effectively detected. 

5 

The system and/or method can analyse for Brugada Syndrome by analysing at least the 
following parameters: PR-duration, ST elevation, ST morphology and T wave mor- 
phology. Hereby, it is achieved that Brugada Syndrome is detected in a very effective 
and fast way. 

10 

The system and/or method can analyse for Right bundle branch block RBBB by ana- 
lysing at least the following parameters: QRS duration, QRS morphology, T wave 
morphology and ST elevation. In this way, an indication of Right bundle branch block 
RBBB is achieved in an effective way. 

15 

The system and/or method can analyse for Left bundle branch block LBBB by analys- 
ing at least the following parameters: QRS duration, R wave morphology and T wave 
morphology. An effective indication of Left bundle branch block LBBB is also possi- 
ble with this method or by using the system. 

20 

The system and/or method can analyse for Short QT Syndrome by analysing at least 
the following parameters: Q-T duration and T wave morphology. Short QT Syndrome 
can also be analysed effectively. 

25 The system and/or method can analyse for Hyperkalemia by analysing at least the fol- 
lowing parameters: P wave morphology, T wave morphology, QRS duration, QT dura- 
tion and PR duration. Hereby, indication of Hyperkalemia is achieved in a highly ef- 
fective way. 

30 The system and/or method can analyse for Hypokalemia by analysing at least the fol- 
lowing parameters: QT duration, T wave morphology and ST depression. Hereby, it is 
achieved that Hypokalemia is indicated effectively. 
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The system and/or method can analyse for peri/myocarditis by analysing at least the 
following parameters: ST elevation, ST morphology, Q wave morphology and PR 
depression. Hereby, it is achieved that peri/myocarditis is detected effectively. 

5 The system and/or method can analyse for Right Ventricular Hypertrophy by analysing 
at least the following parameters: Q wave morphology, QRS duration, S wave mor- 
phology and T wave morphology. Herby, it is achieved that Right Ventricular Hyper- 
trophy is indicated effectively. 

10 The system and/or method can analyse for Left Ventricular Hypertrophy by analysing 
at least the following parameters: Q wave morphology, QRS duration, S wave mor- 
phology and T wave morphology. Also Left Venticular Hypertrophy can be detected 
effectively. 

15 The system and/or method can analyse for Arrhythmogenic Right Ventricular Dyspla- 
sia by analysing at least the following parameters: QRS duration, S wave morphology 
and T wave morphology. Hereby, it is achieved that Arrhythmogenic Right Ventricu- 
lar Dysplasia is detected in an effective way. 

20 Below is described one possible method and a system to illustrate the invention. 

* 

Abstract 

The Long QT Syndrome is a genetic disorder characterized by abnormal cardiac repo- 
larisation resulting in prolonged QT duration, syncopal episodes and increased risk of 
sudden cardiac death. Mutations in the KvLQTl- and HERG genes account for more 

25 than 90% of all LQTS patients. The QT interval duration is the only ECG-based quan- 
tifier of LQTS used in clinical practice today. However duration is only a gross esti- 
mate of repolarisation and does not allow perfect discrimination between KvLQTl, 
HERG and normal subjects. Studies have shown that T-wave morphology parameters 
are useful discriminators in LQTS, but no single parameter has proven to be sufficient. 

30 In this study we present a novel multivariate discrimination method based on a combi- 
nation of T-wave symmetry-, flatness- and duration parameters. 16 subjects were in- 
cluded in the study - 8 normal, 5 HERG and 3 KvLQTl patients. Genotypes were 
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known for all LQTS patients, but one. Standard 12 - lead ECG's were recorded on 
each subject. An automatic ECG event detection algorithm was implemented. The 
signal was highpass filtered and normalized with respect to the isoelectric level to en- 
sure a stable baseline. 4 parameters describing the duration of repolarisation, 6 sym- 
metry- and 1 5 flatness parameters were calculated to characterize each of the T-waves. 
The mean values of lead V5 and the interlead standard deviations were used as pa- 
rameter values. Stepwise discriminant analysis was performed to obtain two discrimi- 

* 

nant functions based on the five strongest discriminatory parameters. The resulting 
discriminant functions include 2 duration-, 2 symmetry- and 1 flatness parameter. The 
two functions classify all subjects correctly (p> 0.0001, p<0.005). Further discriminant 
analysis with a reduced number of parameter categories implied that superior classifi- 
cation is obtained when using all three parameter categories presented. A combination 
of parameters from the three categories symmetry, flatness and duration of repolarisa- 
tion was sufficient to correctly classify ECG recordings from the KvLQTl, HERG and 
normal subjects in this study. This multivariate approach may prove to be a powerful 
clinical tool. 

1. Introduction 

The Long QT Syndrome (LQTS) represents a hereditary genetic disorder characterized 
by the presence of prolonged QT duration on the ECG, syncopal episodes due to po- 
lymorphic ventricular tachycardia (torsade de pointes), and arrythmogenic sudden car- 
diac death. . 

Mutations involving 6 different genes have been identified in LQTS subjects. These 
mutations result in structural and functional changes in ion-channel proteins and cur- 
rents. The changes are manifest by QT prolongation and morphological gene-specific 
repolarisation patterns. The most prevalent genes affected in LQTS patients are 
KvLQTl and HERG which account for more than 90% of LQTS genotype patients. 
The current study focuses on carriers of these two genes. Although some attempts 
have been made to develop quantitative measures that link different repolarisation 
abnormalities to specific LQTS related channel-opathities these methods have so far 
failed to provide a solid diagnostic yield. In current practice the duration of the QT 
interval is the only widely accepted quantifier of ventricular repolarisation. Yet, it has 



10 

been recognized that the duration of the QT interval is only a gross estimate of repo- 
larisation since T-wave morphology is also important when characterizing the QT in- 
terval. This is evidenced by the fact that approximately 10% of all mutation carriers 
have a norma! Bazett corrected QTc (<440ms) and 40% of KvLQTl and HERG carri- 
5 ers show QTc values between 410-470 ms that overlap with non-carriers. Conversely 
only 2% of all carriers present with a normal ST-T pattern and a normal QT interval. 
Morphological aberrations thus carry major implications for the identification of ab- 
normal repolarisation and have been included as diagnostic criteria equivalent to that 
of a positive family history for LQTS. 

10 

Studies have shown that affected KvLQTl patients generally show broad based T~ 
waves with a normal to relatively high amplitude and often without a distinct T - wave 
onset. For individuals with mutations involving the HERG gene the aforementioned 
studies have generally found low amplitude T-waves with bifid T-waves in 60% or 
1 5 more of the carriers. 



Cardiologists already include a qualitative assessment of T-wave morphology from the 
ECG in order to obtain information that augments the clinically established QT inter- 
val measurement and facilitates discrimination between LQTS genotypes. However 
20 qualitative description of repolarisation morphology may be biased due to intra- and 
interpersonal variability thus indicating the need for a standardized quantitative meas- 
ure of this parameter. 

In the following is presented a novel multivariate categorization method that allows 
25 discrimination between KvLQTl, HERG and normal individuals based on Twave 
morphology recorded from 12-lead ECG's. Hallmark morphological features of T- 
waves reported in the literature for these three groups served as inspiration for select- 
ing three primary T-wave characteristics to be assessed. These characteristics are 
symmetry, flatness and duration. 
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2. Methods 



2.1 Subjects 

The study included ECG recordings from 8 female and 8 male subjects. The subjects 
were divided into four groups; 3 KvLQTl (aged 20-48, 2 females), 5 HERG (aged 13- 
5 76, 2 females), 8 normal (aged 23-31, 4 females). Genotypes were known for all 
KvLQTl and HERG subjects with a single exception: 1 patient was categorized as a 
KvLQTl subject by anamnesis and ECG-analysis. In the normal group there were no 
reports of prior cardiac diseases or LQTS family precedent. 

2.2 Data collection 

10 Data acquisition was carried out with the subjects resting in supine position. The 
equipment used for data acquisition was a portable digital ECG recording system, 
"Cardio Perfect Resting ECG system" manufactured by Cardiocontrol. Recording was 
divided into three sessions. Data was collected from 8 leads (I-m, V2-V6) with a 
sampling rate of 1200Hz. Signal recording length was 75 s. in the first session and 150 

15 s. in the last two sessions. 

Following data acquisition, SCP files generated by the Cardio perfect software were 
exported from a MSDE/SQL7 server and subsequently converted to .MAT files using 
SCP-Batch Converter. 

20 

2.3 Algorithm for detection of events in the ECG 

To facilitate evaluation of the repolarisation process and the QT interval, several 
events in the ECG were detected (Qstart, Rtop, Tstart, Ttop and Tend). An algorithm 
for.detecting these events was implemented in Matlab 6.0. 

25 

The method is based on prior work published by Laguna et al. and uses adaptive 
thresholding techniques applied to a digitally filtered and differentiated signal A mi- 
nor extension to the algorithm was incorporated to enable the detection of Tstart. 
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Tstart was detected with a technique equivalent to the technique for detecting Tend. 
Figure 1 shows an example of the result of the event detection algorithm. 

Figure 1. Important events that are used to describe repolarisation are marked by dots 
5 by the event detection algorithm. The algorithm is able to detect the events on all 8 
recorded leads. 

2. 4 Preliminary signal processing 

Evaluation of the QT interval and the repolarisation process was done on the basis of 
an ECG signal with stabilized baseline. This was achieved through preliminary signal 

10 processing. The "raw" ECG was filtered by a minimum order equiripple high pass 
filter with a cut-off frequency of 2 Hz, 60 dB damping in 0 Hz and I dB ripple in the 
passband (N-1276). After filtering, the signal had an almost stable baseline. In order 
to improve stability, isoelectric lines in the signal were estimated from one P-Q inter- 
val (Qstart minus 20 ms) to the following P-Q interval (Qstart minus 20 ms). The sig- 

15 nal was then normalized by subtracting the line value from the corresponding signal 
values. This process is shown in figure 2. 

2. 5 T-wave morphology parameters 

In order to characterize the T-wave morphology, a number of parameters were se- 
lected. The parameters were chosen to cover each of the three categories: Twave 
20 symmetry, T-wave flatness and duration. The parameters are listed and described in 
table 1. 

Parameters S1-S4 and F1-F8 is based on the calculation of modified skewness and 
kurtosis measures. Inspired by the summary measures of probability distributions used 
25 in the field of statistics the T-waves were modelled as probability mass distributions 
(figure 3) and assigned a centre (mean), width (standard deviation), an asymmetry 
measure and a convexity measure. Asymmetry and convexity calculations were then 
carried out based on the modified skewness and kurtosis measures (3 rd and 4 th order 
moments) as follows: 

30 

The total area under the signal, mO, was calculated: 



13 



n=0 



The signal was normalized by the value of the area, mO: 
v{n] = v[n]/m 0 

5 

Normalization facilitated the calculation of the moment functions, since w[n] shares a 
fundamental property with the probability mass function: A total area of 1 . 
The 1 st order moment, ml, was calculated, ml is the mean of the signal: 

AM 

10 m l = ^L,n*w[ri] 

■ 

The 2 nd order moment, m2, was calculated. m2 is the standard deviation of the signal: 
m 2 = (X("-w,) 2 *M'W)" 2 

15 

Figure 2 Isoelectric lines (dashed lines) in the signal are calculated from one P-Q in- 
terval to the following P-Q interval (Qstart - 20 ms). The line values are subtracted 
from the corresponding ECG signal values giving the distances v(n). The result of this 
procedure is shown as an area plot with basis on the zero-line. 

20 Figure 3. a) Example probability mass distribution used when calculating standard 
skewness and kurtosis measures, b) Modified frequency distribution used in this study 
for calculating the modified skewness and kurtosis measures. Signal values v(n) are 
shown in figure 2. 

4/8 



25 PARAMETER DESCRIPTION 



Symmetry 

S 1 Skewness evaluated from Tstart to Tend. 



14 

52 Skewness evaluated from Tstart to Tend with Ttop as mean. 

53 Skewness evaluated in a symmetric interval, 10 % of the Tstart-Tend interval 
surrounding Ttop with Ttop as mean. 

54 Skewness evaluated in a symmetric interval, 20 % of the Tstart-Tend interval 
surrounding Ttop wit Ttop as mean. 

*% 

55 Ratio of the time interval from Tstart to Ttop and the time interval from Ttop to 
Tend. 

56 Ratio of the average slope from Tstart to Ttop and from Ttop to Tend. 
Flatness 

Fl Kurtosis evaluated from Tstart to Tend. 

F2 F 1 normalized by the absolute Rtop-Qnadir value. 

F3 Kurtosis evaluated from Tstart to Tend with Ttop as mean. 

F4 F3 normalized by absolute Rtop-Qnadir value. 

F5 Kurtosis evaluated in a symmetric interval, 10 % of the Tstart-Tend interval sur- 
rounding Ttop with Ttop as mean. 

F6 F5 normalized by absolute Rtop-Qnadir value. 

F7 Kurtosis evaluated in a symmetric interval, 20 % of the Tstart-Tend interval sur- 
rounding Ttop with Ttop as mean. 

F8 Kurtosis normalized by the value of Rtop with Ttop as mean. 

F9 Ratio of the total area under the T-wave from Tstart to Ttop and the correspond- 
ing time interval. 

F 1 0 F9 normal ized by absolute Rtop-Qnadir value. 

Fl 1 Ratio of the total area under the T-wave from Ttop to Tend and the correspond- 
ing time interval. 

F12 Fl 1 normalized by absolute Rtop-Qnadir value. 

F13 Ratio of the total area under the T-wave from Tstart to Tend and the correspond- 
ing time interval. 

F 1 4 F 1 3 normalized by absolute Rtop-Qnadir value. 
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Fl 5 Ratio of the height of Rtop and the width of the Tstart-Tend interval 



Duration 

QTc The Q-T interval normalized by the square root of the R-R interval according to 
Bazett's formula. 

5 D2 Time interval from Tstart to Tend. 
D3 Time interval from Tstart to Ttop. 
D4 Time interval from Ttop to Tend. 

The table above shows a Complete list of the parameters used to characterize T-wave 
10 morphology. Parameters belong to one of three categories: symmetry, flatness and 
duration. 

The 3 rd order moment, m3, was calculated. m3 is the modified skewness of the signal: 

15 

/ « 3 =(|;( M -/» 1 ) 3 *w[»]) ,/3 



Finally the 4 th order moment, m4, was calculated. m4 is the modified kurtosis of the 
20 signal: 



2.6 Data analysis in Matlab 

25 The T-wave morphology parameters for the acquired, pre-processed ECG recordings 
were evaluated using Matlab 6.0. Only valid data were analyzed - i.e. data from leads 
where the signal was not corrupted by high frequency noise and where the event detec- 
tion algorithm was successful in detecting the relevant events with satisfactory preci- 
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sion. Parameter means and standard deviations were calculated for every T-wave in 
the signal on all leads. A great interlead variation in T-wave morphology may be an 
indicator of LQTS. Interlead variance was therefore examined by calculating the stan- 
dard deviation of the lead means for each parameter. 

5 

Only the parameter means from lead V5 and interlead standard deviations were used 
as final parameter values. Hence, for every parameter in table 1, two parameters were 
calculated - one with index "meanVS" and one with index "std" e.g. FlmeanV5 and 
Flstd. 

10 2.7 Statistical analysis 

In order to characterize and classify data from the three groups (KvLQTl, HERG and 
normal), the evaluated parameter values were processed using discriminant analysis. 
The analysis was carried out in SPSS version 11.5. The objective of the discriminant 
analysis was twofold: finding parameters that most efficiently discriminate between 
15 the groups and reducing the number of variables. Therefore a stepwise procedure was 
used with the Mahalanobis D2 as the most appropriate distance measure. 

The entry/removal-criteria were adjusted in order to reduce the number of variables in 
the discriminant functions to achieve a 1:3 ratio between the number of variables and 
20 the population size (N=16). The criteria were empirically chosen to be pentry = 0.045 
and premoval = 0.09 providing the desired 5 variables in the discriminant functions. 

m 

Figure 4. Scatterplot showing classification of individuals by genotype. Separation of 
groups was carried out by 2 discriminant functions with 5 variables that characterize 
25 repolarisation by computation of symmetry, flatness and duration. 

3. Results 

« 

The discriminant functions were based on data from all KvLQTl, HERG and normal 
subjects. The 5 parameters included in both discriminant functions are listed in table 2. 

The discriminative efficiency of both generated functions was statistically significant 
30 after inclusion of all 5 parameters (function 1: p<0.0001, function 2: p<0.005). 
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Variables Entered 



Step 


Entered 


1 


Fllstd 


2 


QTcmeanV5 


3 


S5meanV5 


4 


D4std 


5 


S4meanV5 



5 Table 2. Variables used by the two discriminating functions. Stepwise introduction of 
more variables improved the ability of the functions to discriminate between KvLQTl, 
HERG and normal. 

A scatterplot was generated from the discrimination fimctions and groupings of indi- 
10 vidual genotypes can be seen in figure 4. The dotted lines were read from the SPSS 
generated territorial map and manually added. The lines reflect borderlines where the 
differences between each pair of discrimination functions are zero. All 16 processed 
ECG's were correctly classified and showed at least one discriminatory characteristic 
as defined by the 5 parameters included in the discrimination functions. Cross valida- 
1 5 tion of both discriminant functions was done with the leave-one-out method and all 1 6 
subjects were again correctly grouped. Reducing the number of variables resulted in 
misclassified cases due to lack of one or more discriminatory characteristics. In light 
of this finding we elected to perform further analysis of the selected parameters in or- 
der to investigate the individual contributions of each variable to the separation of the 
20 three primary groups of subjects. Extreme values for all parameters were identified 
and the mean was computed. 

The result is plotted in figure 5. As expected the extent of interlead flatness variation 
observed in HERG and normal individuals was lower than that found in KvLQTl sub- 
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jects. This is evidenced by the Fllstd parameter in figure 5a. When evaluating pa- 
rameter values S4meanV5 and S5meanV5 (figure 5b, d) the extent of 6/8 asymmetry 
in KvLQTl and normal was generally less than that of HERG individuals. Both 
S4meanV5 and S5meanV5 are symmetry parameters and asymmetry in HERG indi- 

5 viduals was augmented in two ways: When bifid T-waves were present the interval 
from Tstart to Ttop was prolonged due to the definition of Ttop used in this study (the 
last highest point on the T-wave). Also, when the initial portion before Ttop was pro- 
longed in HERG individuals better discrimination was possible. Both phenomena 
were observed in HERG subjects. Generally the Bazett corrected QTc observed in 

10 HERG and KvLQTl was higher than that of normal individuals (figure 5e). However 
overlap existed between all three groups preventing separation of the groups by QTc. 
Since no single parameter included in the discrimination functions was able to sepa- 
rate KvLQTl, HERG and normal, we proceeded to investigate the classification effi- 
ciency provided by the three primary categories, represented by the parameters in the 

15 functions. This was carried out by generating new discrimination functions using pa- 
rameters from one category only while excluding the other two. Then, from the new 
discrimination functions three additional functions were generated, this time allowing 
the inclusion of parameters from combinations of two categories. Scatterplots illustrat- 
ing the results of this analysis are shown in figures 6a-f. The first two functions (figure 

20 6a) included parameters that characterize the symmetrical properties of the Twave. 
83.1% of the 16 subjects were correctly classified. Arrows in figure 6a indicate the 3 
misclassified subjects. A second discriminant analysis was performed using flatness 
parameters. This resulted in 93.8% correctly classified subjects. Only one subject was 
not correctly classified as indicated by the arrow on figure 6b. The misclassified case 

25 was the same HERG subject incorrectly classified using symmetry parameters. The 
discriminatory efficiency of duration parameters was also evaluated. Discrimination 
analysis resulted in 93.8% correctly classified subjects. One HERG subject was mis- 
classified as KvLQTl. QTc was 416ms and the subject showed relatively peaked T- 
waves similar to those found in KvLQTl. However the duration parameters failed to 

30 identify this morphological feature, thus reducing classification performance. 

It can be noted that improved classification was obtained using flatness or duration 
parameters versus symmetry parameters and it seemed reasonable to investigate if fur- 
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ther classification improvement could be achieved using a combination of several pa- 
rameter categories. 

Figure 5. a) Fllstd -Interlead standard deviation of the ratio between the total area 
5 under the T-wave from Ttop to Tend and the corresponding time interval, b) 
S4meanV5 - Lead V5 mean modified skewness evaluated in a symmetrical interval 
surrounding Ttop and corresponding to 20% of the interval between Tstart-Tend. c) 
D4std - Interlead standard deviation of the time interval from Ttop to Tend, d) 
S5meanV5 - Lead V5 mean of the ratio between the time interval from Tstart to Ttop 
1 0 and the corresponding time interval from Ttop to Tend, e) Lead V5 mean QTc. 

Figures 6d«f show the results of three separate discriminant analysis using combina- 
tions of parameters from two categories. It can be noted that classification of subjects 
was perfect in all cases, even when repolarisation duration was not considered (figure 
15 6d). 

4. Conclusion and discussion 

The initial discriminant analysis performed in this study resulted in perfect classifica- 
tion of all KvLQTl, HERG and normal subjects. In table 2 it was noted that the dis- 
criminant functions included parameters from all three categories; T-wave symmetry, 
20 T-wave flatness and duration. This is in agreement with the initial hypothesis that a 
combination of repolarisation duration and T-wave morphology characteristics could 
improve discrimination between KvLQTl, HERG and normal. 

To understand why some subjects were misclassified using a reduced set of parameter 
categories (figures 6a-c) the duration parameters and morphological characteristics of 
25 all 16 ECG's were examined. 

Using only symmetry parameters, 3 subjects were misclassified. However no obvious 
visual characteristics on the three misclassified ECG's could be identified that ex- 
plained the incorrect classifications. The Bazett corrected QTc was 347ms for the 
normal subject, 425ms KvLQTl, 476ms HERG. Although an obviously prolonged 
30 QTc was present in the misclassified HERG subject it was not identified using sym- 
metry parameters alone. 
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Discriminant analysis using parameters from the flatness category resulted in only 1 
reclassification. Again no visual characteristics were identified to account for the 
misclassification. Although it was anticipated that the 

Figure 6. a) The result of discriminant analysis using symmetry parameters resulted in 
three misclassified cases (arrows). Visual inspection of the ECG's revealed no appar- 
ent abnormalities to indicate the reason for incorrect misclassification. b) The result of 
discriminant analysis using flatness parameters. One incorrectly classified HERG sub- 
ject was identified (arrow) even though no obvious visual abnormality indicated a dif- 
ferent genotype, c) Result of discriminant analysis using duration parameters. This 
result illustrates the failure of duration parameters to discriminate between KvLQTl, 
HERG and normal (arrow), d-e) Combinations of parameters from two categories il- 
lustrate the improvement in classification efficiency when compared to figures 6a-c 
evaluation of T-wave flatness would be able to discriminate HERG from KvLQTl 
subjects this was not accomplished by using flatness as a single descriptor of repolari- 
sation. Performing discriminant analysis based on the QTc parameter as the only vari- 
able resulted in 1 misclassification. This was not unexpected since it is well known 
that a substantial overlap in QTc values can exist between normal and affected indi- 
viduals. The lack of unambiguous discrimination between all groups by use of the 
QTc parameter alone emphasizes the hypothesis that additional parameters are needed 
to classify LQTS individuals. By combining parameters from two categories it was 
found that the discriminatory strength was increased. 

(figures 6d-f) This was evidenced by the fact that no subjects were misclassified using 
two categories. A particularly interesting finding, was the perfect separation of all sub- 
jects that was obtained using symmetry and flatness parameters with no duration pa- 
rameters included. This result implies the discriminatory strength inherent in parame- 
ters from those two categories. In addition it was found that symmetry or flatness pa- 
rameters combined with duration parameters yielded perfect discrimination between 
all groups. Results from the discriminant analysis using one and two categories indi- 
cate that a combination of more parameter categories strengthen the overall discrimi- 

* 

natory power of the classification functions. Combining these findings with the results 
from the three category discriminant analysis initially performed, it is reasonable to 
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speculate that a substantially improved discrimination between KvLQTl, HERG and 
normal is possible using all three categories of parameters. 

In light of the results obtained in this study we propose a new technique for discrimi- 
5 nating between KvLQTl, HERG and normal subjects. Through multivariate discrimi- 
nant analysis it was found that a combination of two duration parameters and three T- 
wave symmetry-and flatness parameters was sufficient to classify each of the 16 study 
subjects into one of the three distinct groups. Although no single parameter had the 
necessary discriminatory strength to classify the subjects, the combination of multiple 
10 parameters in two discrimination functions was statistically significant (function 1: 
pO.OOOl, function 2: p<0.005). The encouraging results of multivariate repolarisation 
analysis found in this study support the use of symmetry-, flatness- and duration pa- 

» 

rameters to classify LQTS patients. 

15 Further effort must be made to strengthen the statistical impact of this study and to 
investigate whether the proposed symmetry-, flatness- and duration parameters are the 
best suitable for the discrimination functions. 

The use of the proposed multiple parameter categories to classify KvLQTl and HERG 
20 genotypes may prove to be a powerful clinical tool in the making. 

The invention will in the following be described in detail with reference to the draw- 
ing, where figures 1 to 6 are mentioned in the text, and they are as such not further 
described. 

25 

Fig. 7 shows examples of ECG-curvature having ST-elevation myocardial infarction, 
where the active parameters for indicating ST-elevation myocardial infarction is ST- 
elevation, ST-morphology, T-wave morphology and Q-wave morphology. 

30 Fig. 8 shows curvaturic non ST-elevation myocardial infarction, where parameters of 
ST depression, T-wave morphology or Q-wave morphology could be used. 
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Fig. 9 shows a ECG-curvature indicating Cardiomyopathia where the following pa- 
rameters can be used for indication: P-wave morphology, QRS duration, S-wave mor- 

♦ 

phology and T-wave morphology. 

Fig. 10 shows in the same way indications for Brugada Syndrome. Here, the effective 
parameters for indication are PR-duration, ST-elevation, ST morphology and T-wave 
morphology. 

Fig. 1 1 shows curvatures referring to Right bundle branch block RBBB where the pa- 
rameters could be QRS duration, R-wave morphology, T-wave morphology and ST- 
elevation. 

Fig. 12 shows curvatures indicating Left bundle branch block LBBB where the effec- 
tive parameters could be QRS duration, R-wave morphology and T-wave morphology. 

Fig. 13 shows curvatures indicating Short QT syndrome where parameters for indica- 
tion could be Q-T duration and T-wave morphology. 

Fig. 14 shows curvatures indicating Hyperkalemia where parameters effective for in- 
dications are P-wave morphology, T-wave morphology, QRS duration, QT duration 
and PR duration. 

Fig. 15 shows curvatures indicating Hypokalemia where the effective parameters for 
indication seem to be QT duration, T-wave morphology and ST depression. 

Fig 16 shows curvatures indicating Pericarditis where the effective parameters are ST 
elevation, ST morphology, Q-wave morphology and PR depression. 

Fig. 17 shows a curvature indicating Right Ventricular Hypertrophy (RVH) where the 
effective parameters for indication are Q-wave morphology, QRS duration, S-wave 
morphology and T-wave morphology. 
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Fig. 18 shows curvatures indicating Left Ventricular Hypertrophy (LVH) where the 
effective parameters for indication are Q-wave morphology, QRS duration, S-wave 
morphology and T-wave morphology, 

5 Fig. 19 shows curvatures indicating Arrhythmogenic Right Ventricular Dysplasia 
where the parameters are QRS duration, S-wave morphology and T-wave morphology 
is used. 
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ABSTRACT 

The present invention relates to a system or a method for analysing ECG curvature 
where at least one among a number of different parameters is isolated, which system 
has a input means connected to an ECG source, where the different parameters of a 
5 received ECG curvature are indicated and/or isolated and for indicating possible 
symptoms which relate to or are indications of certain deceases, where said deceases 
are known to influence the ECG curvature. 

The aim of the invention is to achieve a system and a method for objective, fast and 
10 effective indication of a number of symptoms derivable from an ECG curve which 
may be indicative of one or more diseases. 

This can be achieved with the system or method previously described, if a first number 
of selected parameters are combined in at least a first mathematical analysis, where the 

1 5 result of the analysis can be represented as a point in a coordinate system comprising 
at least two axes where the system can compare the actual placement in the coordinate 
system with a number of reference parameters stored in the system for indicating dis- 
eases having influence on the ECG curvature. Hereby, it is achieved that any symptom 
of a disease having an indication (influence) in the ECG curvature can be detected in 

20 an objective, automated and very fast way. The system might be used under field con- 
ditions such as in ambulances or in other situations where a fast indication of heart 
diseases is needed in order to help the patient in a correct way as early as possible. 

Fig. 4 
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Probability mass distribution 



T-wave modelled as a probability mass distribution 




•10 12 

Value v(n) 



4 




Q 20 

Tslart 



60 BO 100 120 

sample number 



140 1ED 

Tend 



Fig. 4 



CM 
C 

o 

•s 



7 - 



5 - 



3 - 



1 - 



-1 - 



-3 - 



-5 - 



-7 
-9 




Groups 

O Group Centroids 

<*• HERG 

■ KvLQTI 

° Normal 



Function 1 



Fig. 5a 



,00014 



,00012 



,00010 



tJ 



,00008 



,00006 



[] 



,00004 



,00002 



0,00000 



1 —i 

Normal KvLQTI 



HIGH 
LOW 
D MEAN 



HERG 



Coov orovided by USPTO from the IFW Image Database on 11/18/2004 



Fig. 5b 




Normal KvLQTI HERG 
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ST-elevation myocardial infarction 




• ST elevation 



• ST morphology 

• T-wave morphology 

• Q-wave morphology 

Traditionally the presence of an ST elevation myocardial infarction (STEMI) is in 
the ECG diagnosed in the early phase by the ST elevation in the J-point and the 
subjective evaluation of the shape of the ST-morphology together with the 
progressive development of Q-waves. With the present method the diagnosis will 
be based on the combination of the objective measures and new parameters making 
the subjective changes into objective parameters. E.g the shape of the ST 
segment (ST morphology) will be characterised by I or more parameters, by which 
it is expected that it will be possible to improve the diagnosis in an automated way. 
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Non ST-elevation myocardial infarction 
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• Q-wave morphology 
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Fig. 9 
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Fig. 16 

Pericarditis 
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Fig. 19 



Arrhythmogenic Right Ventricular Dysplasia 
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